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Abstract 

Based on the new type of random walk process called the Potentials of Unbalanced Complex 
Kinetics (PUCK) model, we theoretically show that the price diffusion in large scales is amplified 
2(2 + 6) _1 times, where b is the coefficient of quadratic term of the potential. In short time scales 
the price diffusion depends on the size M of the super moving average. Both numerical simulations 
and real data analysis of Yen-Dollar rates are consistent with theoretical analysis. 
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I. INTRODUCTION 



Crashes and uncontrollable hikes can often occur in financial markets. Such changes 
of the prices confuse the market and damage the economy because they start abruptly in 
many cases. Therefore, techniques to measure the probabilistic risk of sudden change in the 
prices have been studied using tick- by-tick data [1]. Recently, it was empirically found that 
change of prices can be approximated by the Fokker-Planck equation and the new type of 
random walk in a potential field [2,3,4,5]. The potential field is approximated by a quadratic 
function with its center given by the moving average of past market prices. This random 
walk model is called the Potentials of Unbalanced Complex Kinetics (PUCK) model in which 
the potential slowly changes in the market [3,4]. In this paper, we focus on the diffusion 
properties of this random walk process and calculate the diffusion coefficient which is helpful 
for estimating the market risk. 

We first review an empirical derivation of the PUCK model. We next show that the 
statistically steady condition of price fluctuations depends on the potential field, and clarify 
relationships between the price diffusion and the potential field. We finally demonstrate 
that the price diffusion in short time scales depends on the size of moving average, however, 
large scale diffusion properties are independent of the moving average. In the paper, we 
used all the Bid record (about 20 million ticks) of the exchange rates for Yen/Dollar that 
were traded by the term from 1989 to 2002 to find the firm statistical laws. 

II. EMPIRICAL DERIVATION OF PUCK MODEL 

Prices in financial markets always have violent fluctuation in a short time scale. We first 
eliminate the uncorrelated noise rj(t) from the price P(t) in order to reduce the statistical 
error. We next investigate the dynamics of the price. 

We can perform this noise elimination process by introducing an optimum moving average 



P(t): 



P(t) = P(t) + V (t), 



(1) 




(2) 
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where P(t) is a price, rj(t) is an uncorrelated noise and uik gives the weight factors where the 
time is measured by ticks. The weight factors are calculated by using Yule- Walker equation 
[6,7,8]. In a case of Yen-Dollar rate, the weight factor uik follows an exponential function 
whose characteristic decay time is about 30 seconds [6]. 

We investigate a behavior of the optimum moving average P(t) obtained by eliminating 
the uncorrelated noise rj(t) from the price P(t). We introduce a super moving average -Pm(^) 
defined by 

= (3) 

In financial markets, it is found that the following relationship between P(t) and Pm{^) 
holds for a certain range of M[2], 

WTT)-W) = -\- 1 ^-^(W)-P^y + f(t), (4) 

where the mean of noise f(t) is zero. The Eq.(@J) means that the price change can be 
approximated by a random walk in a quadratic potential field whose center is given by the 
moving average of past prices. It is known that the potential coefficient b(t) has a long 
autocorrelation [3,4]. 



III. STATISTICALLY STEADY CONDITION OF PRICE FLUCTUATIONS IN 
THE PUCK MODEL 

We focus on Eq.(JJJ) with the case of a constant b because the coefficient b(t) is known to 
change slowly in financial markets. Eq.pjl is transformed as follows, 



b I 2 \ 

W+T)-W) = ~2 I M(M _ 1} E( M ' fc ) (P(t-k+l)-W^kj)j +f(t). (5) 

This is a type of AR process for price difference when b is a constant. We can estimate the 
conditions of b to make the AR process being statistically steady. Eq.(j3J) is transformed by 
the following determinant, 
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where 



X t = AXt-x + F t 

= (A) 1 Xq + [F t + AF t -i + 



■ + {At 1 F x ) , 



(6) 
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If lim (Af = 0, the time series of P(t + 1) — P(t) is in a weakly steady state because X t 
becomes independent of the initial value of X . This condition is fulfilled when the absolute 
values of all eigenvalues of A are less than 1. For example, when M = 2 the time series 



of P(t+ 1) — P(t) is statistically steady if the potential coefficient satisfies \b\ < 2. When 
M = 3 the eigenvalues of A are given by 



Ai 



«i + \Ja\ + 4a 2 



A, 



ax — y/a{ + 4a 2 



(8) 
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Solving the steady state condition we find that the potential coefficient should be in the 

range of —2 < b < 6 for M — 3. Numerically calculating the eigenvalues for 2 < M < 20 



as shown in Fig.l we find that the time series of P(t + 1) — P(t) is statistically stationary if 
the potential coefficient is in the following range, 



(9) 



-2 < b < 2(M — 1) when M is even number 
—2 < b < 2M when M is odd number 

Outside the condition of Eq.Q, the range of price fluctuations increases indefinitely de 
pending on the time t. 



IV. DIFFUSION OF PRICES IN MARKET POTENTIAL FIELD 



As the potential coefficient b(t) has a long autocorrelation, we can calculate the future 
price diffusion using Eq.pjl. This prediction is crucial in order to evaluate the risks of market. 
We clarify statistical laws of price diffusion described by Eq.(JU) using both simulations and 
theoretical analysis. 

By simulating Eq.(jlJ) for the case of b(t) is a constant, we investigate the standard devi- 
ation on a time scale T defined by 



<T b (T) = J((P(t + T)-P(t) . (10) 



In Fig.2 we plot a b (T) for M = 4, 16, 64, 256 when b = -1.5 and 6 = 2. Here, f(t) is 
the Gaussian random number whose standard deviation is 1. The time scale where the 
Hurst exponent converges 0.5 depends on M. For example, the time scale is about 10 3 when 
M = 256, while the time scale is around T = 20 when M = 4 

We can estimate the standard deviation <Jb{T) in the long time scale limit. The variance 



of change of optimum moving average price P(t) is given by 
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By introducing a rough approximation 

((PTt + Tj -¥({}) ■ (P(t + T-*)-P(R)) « ((^(< + T) -PW) 2 ^, (12) 
we have the following simple formulation after some calculation, 

<T) = L^-\ a b=0 (T), (13) 
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where T 3> M, <Jb(T) = J ^ (p(t + t) - p{t)} j and (Jb=o(T) is the standard deviation when 
6 = 0. In this long time scales, the standard deviation of the price is amplified 2(2 + b)~ x 
times by the potential field. In Fig. 3 we show relationships between the potential coefficient 
b and the ratio of standard deviation cr 6 (T) / ab = o(T) when T = 10 5 by simulating Eq.Q 
for M = 2, 4, 8, 16, 32. We can confirm that the price diffusion of numerical simulations 
follows Eq. (jl3j) independent of M in the long time scale. From Eq. (jl3|) we can theoretically 
find that in the long time scales the price diffusion is independent of M, independent of b, 
a b (T) oc T - 5 . 

V. DIFFUSION OF YEN-DOLLAR RATES 

From dataset of real market prices, we can estimate the value of b, and presume the best 
value of M by comparing the price diffusion of numerical simulations to real price diffusion. 
Fig. 4 shows the diffusions of Yen-Dollar rates from 3:35 to 8:35 and from 9:25 to 23:25 
in 11/9/2001, the day of terrorism. The rates were stable till 8:35 and it became quite 
unstable after 9:25. The rates until 8:35 follow a slow diffusion in short time scales, namely, 
the market has an attractive potential b > 0. As shown in Fig.4, we can recreate the rate 
diffusion in all time scales by simulating Eq.(0J) with b = 0.9 and M = 14 ticks. We next 
focus on the unstable rates after 9:25 that follow a fast diffusion in short time scales. When 
b = —0.8 and M = 14 ticks, the price diffusion of numerical simulation is also consistent 
with the rate diffusion. The price diffusion depends on M in the short time scales, although 
the price diffusion is independents of M in the long time scales. Therefore, if M ^ 14 ticks, 
there are gaps between the price diffusion of numerical simulations and the rate diffusion in 
the short time scales as shown in Fig.4. In other real markets, we can also estimate the best 
value of M from such characteristics of price diffusion. 

VI. DISCUSSION 

We can approximate the change of market prices by the random walk in a potential field. 
The potential field is well approximated by a quadratic function with its center given by the 
moving average of past prices. The random walk process is called the PUCK model. By 
analyzing the model, we clarified that the statistically steady condition of price fluctuations 
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depends on the potential coefficient 5, and we also theoretically proved that the price diffusion 
in the long time scales is amplified 2(2 + b)^ 1 times, independent of the size M of super 
moving average. In short time scales the price diffusion depends on M. We can estimate the 
best value of M in real financial markets by observing this dependence. We recreated the 
diffusion of Yen-Dollar rates in all time scales by the PUCK model. The potential coefficient 
b is helpful to measure the probabilistic risk of sudden change in the prices. We may be able 
to build better financial options that offset the risk by applying the price diffusion of the 
PUCK model. 

Acknowledgments 

This work is partly supported by Research Fellowships of the Japan Society for the 
Promotion of Science for Young Scientists (T.M.). The authors appreciate H. Moriya of 
Oxford Financial Education Co Ltd. for providing the tick data. 



[1] J. P. Bouchaud and M. Potters, Theory of financial risks, Cambridge university press, 2000. 
[2] M. Takayasu,T. Mizuno, T. Ohnishi, and H. Takayasu, Temporal Characteristics of Moving 
Average of Foreign Exchange Markets, in Practical Fruits of Econophysics (Springer Verlag, 
Tokyo. 2006) 29-32. 

[3] M. Takayasu, T. Mizuno and H. Takayasu, Potentials of Unbalanced Complex Kinetics Ob- 
served in Market Time Series, physics/0509020 

[4] M. Takayasu, T. Mizuno and H. Takayasu, Potential Force Observed in Market Dynamics, 
Physica A, to appear. 

[5] R. Friedrich, How to Quantify Deterministic and Random Influences on the Statistics of the 
Foreign Exchange Market, PRL 84, 5224-5227, 2000. 

[6] T. Ohnishi, T. Mizuno, K. Aihara, M. Takayasu, H. Takayasu, Statistical properties of the 
moving average price in dollar-yen exchange rates, Physica A 344, pp. 207-210, 2004. 
[7] T. Mizuno, M. Takayasu, and H. Takayasu, Modeling a foreign exchange rate using moving 
average of Yen-Dollar market data, in Practical Fruits of Econophysics (Springer Verlag, Tokyo. 
2006) 57-61. 



7 



[8] T. Ohnishi, T. Mizuno, K. Aihara, M. Takayasu, H. Takayasu, Systematic tuning of op- 
timal weighted-moving-average of yen-dollar market data, in Practical Fruits of Econophysics 
(Springer Verlag, Tokyo. 2006) 62-66. 



8 



+ M=3 




-1.0 -0.5 0.0 0.5 1.0 



c 1 

FIG. 1: The eigenvalues A of matrix A. The c\ and C2 are real and imaginary parts of A. The 
pluses are A for M = 3 when —2 < b < 6. The squares show A for M = 11 when —2 < b < 22. 
The diamonds express A for M = 20 when — 2 < b < 38. The circle indicates |A| = 1. 




FIG. 2: Standard deviation a(T) of price changes in the time scale T. The heavy lines show 
price diffusions of numerical simulation for M = 4(top), 16(second), 64(third), 256(bottom). Here, 
b = —1.5 in (a) and b = 2 in (b). The standard deviation of f(t) is 1. The guideline indicates 
a(T) oc T - 5 . 
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FIG. 3: Relationship between potential coefficient b and standard deviation ratio <Jb{T) / (Tb=o(T) 
when the time scale T = 10 5 . Theory (dashed line) and numerical simulations for M = 2 (plus), 
4(diamond), 8(square), 16(black triangle), 32(white triangle) . [4] 




T [ticks] 

FIG. 4: Diffusion of Yen-Dollar rates. The squares and triangles show the diffusion from 3:35 to 
8:35 and from 9:25 to 23:25 in 11/9/2001. The top 3 lines are the numerical simulation of Eq.(j2J) 
for M = 2, 14, 100 ticks when b = —0.8. The standard deviation of f(t) is 0.0115 yen/dollar. The 
bottom 3 lines show the numerical simulation for M = 2, 14, 100 ticks when b = 0.9. The standard 
deviation of f(t) is 0.0054 yen/dollar. The straight line expresses the normal diffusion with the 
slope 0.5. 
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